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ABSTRACT 

Distributed massive multiple-input multiple-output (MIMO) 
combines the array gain of coherent MIMO processing with 
the proximity gains of distributed antenna setups. In this pa¬ 
per, we analyze how transceiver hardware impairments affect 
the downlink with maximum ratio transmission. We derive 
closed-form spectral efficiencies expressions and study their 
asymptotic behavior as the number of the antennas increases. 
We prove a scaling law on the hardware quality, which reveals 
that massive MIMO is resilient to additive distortions, while 
multiplicative phase noise is a limiting factor. It is also better 
to have separate oscillators at each antenna than one per BS. 

1. INTRODUCTION 

Cellular radio access networks (RANs) have conventionally 
consisted of one single-antenna base station (BS) per cell that 
served one user equipment (UE) per time-frequency resource. 
Since the increasing data traffic calls for higher spectral effi¬ 
ciencies [bit/symbol/cell], the RAN structure is now evolv¬ 
ing to enable coherent downlink (DL) transmission to multi¬ 
ple UEs per resource symbol. The LTE-A standard has basic 
support for multi-user MIMO using a handful of co-located 
antennas. The massive MIMO concept from [ 1 ] takes multi¬ 
user MIMO to the 5G era by using hundreds of BS antennas 
to serve tens of UEs in parallel on each resource block. 

As the cellular concept is evolving, we can also question 
whether future BSs should be in the cell centers as in the past 
or distributed over the cells. The cloud RAN concept from [ 2] 
provides an efficient way to operate distributed antenna ar¬ 
rays and perform the coherent processing required by massive 
MIMO. While the majority of works on massive MIMO con¬ 
siders co-located arrays, the recent works @0 show that dis¬ 
tributed massive MIMO can provide even higher spectral effi¬ 
ciencies than co-located deployments, due to proximity gains. 

A potential showstopper for distributed massive MIMO 
would be if the technology is too sensitive to transceiver hard¬ 
ware impairments; for example, phase noise in local oscilla¬ 
tors (LOs), amplifier non-linearities, non-ideal analog filters, 
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and finite-precision analog/digital converters. The impact of 
hardware impairments on massive MIMO has received con¬ 
siderable attention in recent years l|5HTT], but only 0 consid¬ 
ered distributed arrays. The paper (6) showed that it is of fun¬ 
damental importance to include hardware impairments in the 
performance analysis, since this can be a main limiting factor 
in systems with many antennas. Nevertheless, @@ showed 
that massive MIMO is resilient to additive distortions origi¬ 
nating from the BS. Multiplicative distortions such as phase 
noise can, however, hinder the system performance. These 
works use simplified stochastic impairment models, but the 
validity of the results has been confirmed in GD by simula¬ 
tions based on sophisticated and realistic models. 

For distributed arrays, an important question is whether 
the antennas should share a common LO (CLO) or if each 
antenna should be equipped with a separate LO (SLO). A 
number of recent works have looked into how this design 
choice impacts the severeness of the phase noise 0HHBD- 
The papers Q|7][8|[T0]] seem to establish the consensus that 
a setup with SLOs is preferable in the uplink (UL), since the 
independent phase rotations average out over the BS anten¬ 
nas. However, the answer is still open when it comes to the 
DL; |9 | showed that a CLO is preferable for non-fading chan¬ 
nels, while |[10) considered fading single-cell systems and 
claimed that CLO prevails for few BS antennas (per user) or 
high SNR, and SLOs are desirable in the opposite cases. 

In this paper, we extend our previous UL work in (5j to 
the DL. We consider a multi-cell massive MIMO system with 
distributed arrays and three kinds of hardware impairments: 
phase noise, distortion noise, and noise amplification. We de¬ 
rive new spectral efficiency expressions for maximum ratio 
transmission (MRT), which establish a performance baseline 
in hardware-impaired multi-cell scenarios. These expressions 
are used to prove how the hardware quality may scale with 
the number of antennas. The analysis shows that SLOs is sys¬ 
tematically a better choice than CLO also in the DL. 

2. SYSTEM MODEL 

We consider a cellular network with L cells that operate in 
a synchronized time-division duplex (TDD) mode. Each cell 



Uplink data 

t G {—tql+Ij • • • ? 0} 


Uplink pilots 

t € {1,... ,B} 


Downlink data 

t G {B+1, • • •, B+tdl} 


Fig. 1 . Illustration of the TDD protocol where each coherence 
block consists of T = tul + tdl + B symbols. 


serves K single-antenna UEs using a BS equipped with N an¬ 
tennas, which can be arbitrarily distributed over the coverage 
area. The TDD protocol divides the time-frequency resources 
into coherence blocks, as illustrated in Fig.[l] Each block con¬ 
sists of T symbols with time indices t = —tul + 1, • • •, B + 
tdl, whereof 7 \jl are UL data symbols, B are UL pilots, and 
tdl are DL data symbols. Note that T = 7ul+Tdl+-B- 
Let (*) T and (*) H denote the transpose and conjugate trans¬ 
pose, respectively. The channel response between UE k in 
cell l and BS j is a constant vector hjik = [hff k ... /i^] T G 

C N within each block, where h^l is the channel response for 
the nth BS antenna. The channels are assumed to be Rayleigh 
fading as 

h jik ~ CM(0, A ji k ), (1) 

where the covariance matrix is Ajik = diag(A..., A^). 

The average channel attenuation A^ > 0 is different for each 
combination of cell indices, UE index, and BS antenna index 
n. It depends, for example, on how the BS antennas are dis¬ 
tributed in the cell and on the UE positions. 

2.1. Uplink Model with Hardware Impairments 

A main goal of this paper is to investigate how transceiver 
hardware impairments impact the DL spectral efficiency. We 
mainly consider impairments at the BSs, since the massive 
MIMO deployment constraints (e.g., cost, power, synchro¬ 
nization, and size restrictions) are likely to lead to BS hard¬ 
ware of lower quality than in contemporary networks. 

To this end, we adopt the UL system model from (5j and 
generalize it to also cover the DL. Since the BSs in massive 
MIMO use channel estimates from the UL to perform trans¬ 
mit precoding in the DL, we need to model both directions of 
the links. As in j5j, the received UL signal y j(t) G C N in 
cell j at symbol time t G { —tul + 1, • • •, B} is modeled as 

L 

yj(t) = . (t) HjiX/ (t) + Vj(t) + rij (t) (2) 

1 = 1 

where x/(t) = [xn (t) ... G C K contains pilot/data 

symbols from UEs in cell l and the channel matrix from these 
UEs to BS j is H ji = [hjn ... h p K \ G C NxK . The sym¬ 
bols from UE k in cell j have power = E{\xjk{t)\ 2 }, 
where E{-} denotes the expected value of a random variable. 

The matrix = diag ..., models 

the multiplicative effect of phase noise (with % = \/—1). The 
variable (j)j n (t) is the phase rotation at the nth BS antenna in 
cell j at time t , and it is modeled as a Wiener process G3 


4>jn(f) ~ — 1), 5) where S > 0 is the variance of the 

phase-noise increments. We consider two implementations: 

1. Common LO (CLO): (t) =... = <f>j n (t) within a cell. 

2. Separate LOs (SLOs): All <pj n (t ) are independent. 

The above represent having one LO that feeds all antennas at 
BS j or one separate LO connected to each of the N antennas. 

Moreover, Vj(t) ~ CA/*(0, Yj (t)) is additive distortion 
noise (e.g., from finite-precision quantization, non-linearities, 
and interference leakage in the frequency domain). It is pro¬ 
portional to the received signal power at the antenna and un¬ 
correlated between antennas (5][l3) : 

TjV) - K uLEZ]^ Ldia s(l /l i;fcl 2 ’---’l /l S ) l 2 ) (3) 

1=1 k=1 ^ ' 

where k ul > 0 is the proportionality coefficient. 

Finally, rjj (t) CJ\f (0, crg S I n) is the receiver noise with 

variance cr| s (including noise amplification in circuits). 

2.2. Downlink Model with Hardware Impairments 

Similar to the UL, we model the received DL signal Zjk{t ) G 
C at UE k in cell j at time t G {B + 1,..., B + tdl} as 

z jk(t) w lm (t)s lm (t) + il> l (t)\+r) jk (t) 

1 = 1 \ 771=1 / 

(4) 

where si m (t ) is the DL data symbol (with power pft = 

E (hmW| 2 }) and W im(t) = K ( m(0 • • • '“’imWY G CN is 

the corresponding linear precoding vector. The receiver noise 
is T]jk(t) CM(0, cr ue)j where cr ^ E is the variance (includ¬ 
ing noise amplification). The phase-noise matrix was 

defined earlier, while ^ (t) ~CM(0, 'f'j) is the additive dis¬ 
tortion in the DL (e.g., due to non-linearities and leakage in 
the frequency domain). Similar to ([3}, the distortion at a cer¬ 
tain antenna is proportional to the transmit power at this an¬ 
tenna and uncorrelated with the distortions at other antennas: 

*3 ~ K DL (kj'fc (*)| 2 . • • • , \Wjk\t) I 2 ) 

k= 1 

where k dl > 0 is the proportionality coefficient. 

This system model is used in the next section to com¬ 
pute achievable DL spectral efficiencies. These depend on 
the level of hardware impairments, as characterized by the 
variance of the phase-noise increments S, the distortion noise 
proportionality coefficients kjjl, kdl, and the receiver noise 
variances cr| s , a^j E . The results are applicable for any 
and for each j and fc, thus under arbitrary power control. 

3. DOWNLINK PERFORMANCE ANALYSIS 

In this section, we derive the DL spectral efficiency per UE 
and study its asymptotic behavior (when N is large) to under¬ 
stand the impact of hardware impairments. 
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3.1. Uplink Channel Estimation 

In order to perform coherent transmit precoding in the DL, 
each BS acquires the channels to its UEs by using the UL 
pilots. The pilot sequence of UE k in cell j is defined as 
Xjfc = [xjk( 1) • • • Xjk(B)] T G C Bxl . The analysis in this 
paper holds for arbitrary pilot sequences (with \xj k (b)\ 2 = 
p B k for b = 1 ,..., B ), while we consider columns from a 
Fourier matrix in Sec. [4] (to achieve mutual orthogonality and 
constant energy per symbol). Since the effective channels 

= E) 0 ^£)h. jik (5) 

depend on the phase-noise and are different at every symbol 
time t , we need a channel estimator that provides new esti¬ 
mates at each t. Such a linear minimum mean-squared error 
(MMSE) estimator was derived in (5) and is described below, 
since the expressions are used in the forthcoming analysis. 
Lemma 1 (Th. 1 in [5]). Let — [yj(l) ••• y ){B)] T G 
C BN denote the combined received signal in cell j from the 
pilot transmission. The LMMSE estimate of h ji k (t) at any 
symbol time t G { —tul + 1, • • •, B + tdl} for any l and k is 

h jik(t) = (x" fc D s(t) <g> A jik) ( 6 ) 

where 0 denotes the Kronecker product, 

Di ( t) =diag^e“ 5 l*- 1 l ) ... )e -fl t - B l^ ) ( 7 ) 


Lemma 2. Suppose that UE k in cell j knows the channel and 
interference statistics, but not the channel realizations. An 
achievable lower bound on the ergodic capacity of this UE is 

-j -B + TDL 

R ik = T l°g 2 (l + SINR, fc (*)) [bit/symbol] (11) 

t=B -\-1 

where SINR j k (t) is given in (T 3 | ) at the top of this page. 

Proof. As in [7), we compute one spectral efficiency for each 
t G {B+ 1, 5+tdl} since the effective channels vary with t. 
The expression is obtained by using the signal received over 
the average channel ¥,{hP-- k (f)ujjk(t)} for decoding, while 
treating the signal received over the uncorrelated deviation 
from this average value, the inter-user interference and distor¬ 
tion noise as worst-case Gaussian noise in the decoder. □ 

The expression in Lemma [2] is a reasonable bound on the 
practical performance that can be achieved using simple sig¬ 
nal processing at the UE (i.e., detect the useful signal and 
treat everything unknown as Gaussian noise). The SINR ex¬ 
pression in ( p~3] ) contains a number of expectations that can be 
computed numerically for any choice of precoding vectors. 
Next, we provide closed-form expressions for MRT. 

Theorem 1. If MRT is used, then the expectations in SINR^ ( t ) 
of Lemma [2] are computed as in (p~4])—(p~5]) at the top of the 
next page (where e n denotes the nth column of I at). 


J 61,62 


L K 

MEE 0 A jUm T cTggl^jv? (^) 

£=1 m= 1 

and the element at position (bi, 62 ) in X^ m G C BxB is 

frf J m( 1 + K UL)> h=b 2 , 

{x im (T hl )x} m (T b2 )e~bi ^ b 2 . 

(9) 

The corresponding error covariance matrix is 

= Aj/fc —(x^D^) 0 (D T s ^-xi k <g> Ajik). 

Proof. This lemma follows from adapting results in 15) to the 
notation and power constraints considered herein. □ 


3.2. Downlink Spectral Efficiency 

Next, we derive achievable DL spectral efficiencies, using 
normalized linear precoding vectors of the general form 

Ujk(t) 


Wjk(t) v^MFT 

With this notation, MRT is given by Wj k (t) = h jj k (t). 


( 10 ) 


Proof. Follows from straightforward computation of the ex¬ 
pectations, whereof some are the same as in [5] Th. 2]. □ 

3.3. Asymptotic Behavior and Scaling Laws 


Next, we investigate the behavior at large N. For tractability, 
we consider A < 00 spatially separated subarrays each with 
^ antennas. Recall that these antennas are either controled 
by a common LO that sends clock signals or separate LOs at 
each antenna. The channel covariance matrices then factorize 
as 

A jik = R jtk ® % ( 16 ) 

where = diag(A^,..., A^) eC AxA and Ais the 
average channel attenuation between subarray a in cell j and 
UE k in cell l. By letting the number of antennas in each 
subarray grow large, we obtain the following property. 
Corollary 1. If MRT is used and the channel covariance ma¬ 
trices can be factorized as in ( p~ 6 ]> , then 


SINR jk (t) 
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(17) 
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if SLOs 
ifCLO 
if SLOs 


where the signal part is 

S jk = tr ((x^Dj W ®AjJ)$ i 1 (D 6(t) x ife 0A^)) 

with = E£=i Em=i x fm ® A^ + crlgl^B, where the 
interference terms Timjk with a CLO are 

E E tow®<) 

jCLO _ Qi = lQ 2 = l _ 

tr 0 A «m) ( D Js(t)*?fc 0 A Uk)) 

X*i 1 ((Xjfe-«uL 4 Ll B)®e 0 l e" 2 )$, 1 (D^ (t) x im 0e a2 ) 

and the interference terms with SLOs are 

= tr ((xL D <5(t) ® A i*m) 0 A/Jfe)) . 


The notation 0(-^) is used for terms that go to zero as or 
faster when N —>■ oo, while e a is the ath column of 1^. 

Proof. Follows from dividing all terms in SINR jk(t) by ^ 
and then analyze the expressions for MRT in (p~4])—(p~5]) . □ 

This corollary does not contain /^ul, ^dl, &bs> or a UE’ 
thus it shows that the impact of distortion noise and receiver 
noise vanishes as N oo. The asymptotic SINRs are only 
limited by the channel distributions, pilot-contaminated inter¬ 
ference, and phase noise. This means that distributed massive 
MIMO systems can handle larger additive distortions than 
conventional systems, as manifested by the next corollary. 

Corollary 2. Suppose that Avttt ^Avttt rj-ZV L ^f>t o-ZV" L 

°BS = °bs,o A22 . O'UE = °ue,o N * 2 > and 5= <5 0 (1 + ln(7V 23 )), 
for some scaling exponents £ 1 , 22,23 > 0 and constants 
^ul,o, ^dl ,05 ^bSjO’ ^UEjO’ ^0 ^ 0- The SINRs, SINR,j&(£), 

with MRT converge to non-zero limits as N oc if 

{ max( 2 i, 22 ) < \ and 23 = 0 for a CLO 
max( 2 !, 2 2 ) + 2 3 6 °\ Td ^- b \ < i for SLOs. 

Proof Based on (p4|)-(p~5]); see (5j for a similar proof. □ 


This corollary shows that the DL can handle additive dis¬ 
tortions with variances that scale as \f~N (i.e., 21 =22 = |), 
while achieving decent performance. The scaling law also 
shows that the phase noise variance with SLOs can increase 
logarithmically with N, while this is not allowed with a CLO. 
This proves that massive MIMO with SLOs are preferable in 
the DL, at least when the number of antennas is large. The 
scaling law holds also for precoders that are better than MRT. 

4. NUMERICAL RESULTS 

The analytic results are corroborated for the distributed mas¬ 
sive MIMO setup in Fig. [2] This is a wrap-around topology 
with 16 cells of 400 x 400 meters, each consisting of A = 4 
subarrays with ^ antennas located 100 meters from the cell 
center. The K = 15 UEs per cell are uniformly distributed, 
with a minimum distance of 25 meters from the subarrays. 
The transmit powers are p® k L = p^ k = — 50 dBm/Hz for all j 
and k (e.g., 100 mW over 10 MHz). The channel attenuations 

are modeled as in [5]: A^ = 10 s ^ fc)_ 1 , 53 /(d^) 3 ' 76 , where 
d^l is the distance in meters between BS antenna n in cell j 

and UE k in cell l and s^ k ~ jV(0,3.16) is shadow-fading. 

The hardware impairments are characterized by the dis¬ 
tortion proportionality coefficients ^ul = ^dl = 0.03, the 
variance of phase noise increments S = 1 • 10 “ 5 , and the re¬ 
ceiver noise powers cr| s = a EE = —169 dBm/Hz (with 5 dB 
noise amplification). These are also the initial constants when 
we scale the hardware quality based on Corollary [2] 

Fig. [3] shows the average spectral efficiency per UE. The 
coherence block contains T = 300 symbols, whereof 5 = 15 
symbols are used for pilot sequences and tdl = 285 for DL 
payload data. Hardware impairments incur a performance 
loss as compared to ideal hardware. The gap is small with 
SLOs, but larger with a CLO. This validates the analytic ob¬ 
servation that SLOs is the better choice in massive MIMO. 

The figure also illustrates the scaling law established by 
Corollary [2] The middle curves show the behavior when sat- 





400 meter 



Fig. 2. Illustration of the multi-cell distributed massive MIMO 
scenario considered in the numerical evaluation. 


isfying the scaling law (z i = Z 2 = 0.48 with a CLO and 
adding also 2:3 = 0.48 with SLOs). By gradually degrading 
the hardware with N, there is a performance loss at every N, 
but the curves are still increasing with N. The performance 
loss is small for SLOs, but very large for a CLO. The curves at 
the bottom are for a case when the scaling law is not satisfied, 
which gives a performance that goes to zero as N —>► 00 . 


5. CONCLUDING REMARKS 

We have analyzed the DL performance of distributed massive 
MIMO systems, with focus on the impact of hardware impair¬ 
ments. We have proved that additive distortions have smaller 
impact on massive MIMO than conventional networks, since 
the variance may increase as with little performance loss. 
Multiplicative phase noise can be more severe, but the perfor¬ 
mance is better if each BS antenna has a separate oscillator. 

The DL analytic results in this paper are in line with pre¬ 
vious UL results in mmm- This is natural since the UL- 
DL duality for systems with linear processing implies that 
the same performance is achievable in both directions (if the 
power allocation is optimized). However, our results stand 
in contrast to the recent works (9]|T0) where the DL behave 
differently than the UL when it comes to phase noise. This 
is due to different system models: © considers high SNRs 
in non-fading single-user cases, while ||T0]| considers a single 
cell with relatively good CSI. In comparison, we consider a 
generalized multi-cell setup with more inter-user interference 
and thus lower SINRs. 
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